function H  = calculateH(n)
m=2*n+1;
H=zeros(2*n-1);
for u=-n/2:n/2-1
    for v=-n/2:n/2-1
        for s=-n/2:n/2-1
            for t=-n/2:n/2-1
                tmp=0;
                for k=-n:n
                    for l=-n/2:n/2    
                            tmp=tmp+exp(2*pi*1i*k*((u-s)-2*l*(v-t)/n)/m)+...
                                exp(2*pi*1i*k*((v-t)-2*l*(u-s)/n)/m);
                    end
                end
                H(u-s+n,v-t+n)=tmp;
            end
        end
    end
end
   